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Abstract 

We calculate ground-state energies and densities of a helium atom confined in an impenetrable 



I 



O spherical box within density functional theory. These calculations are performed by variationally 
^ solving Kohn-Sham equation with the ground-state orbital expanded in terms of Slater-type or- 

• ^ bitals. Using the ground-state densities we then calculate static linear polarizability and nonlinear 
^hyperpolarizability and study their variation with the radius of confinement. We find that polar- 

I— i],zability decreases monotonically with decreasing confinement radius and the hyperpolarizability 
not only decreases but also undergoes a change in sign in the strong confinement regime. 

m 

^PACS numbers: 31.15Bs, 31.15Ew, 36.40Vz, 

o 
o 
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I. INTRODUCTION 



The properties of atoms and molecules undergo drastic change when they are spatially 
confined in either penetrable or impenetrable cavity as compared to their free counterparts. 
In recent years the topic of confined atoms has been attracting lot of attention and it has 
become a field of active research The main reason for this interest in spatially confined 
atoms and molecules is their applicability to several problems of physics and chemistry. For 
example, the model of an atom confined in an impenetrable sphere has been employed to 
simulate the effect of high pressure on the physical properties of atoms, ions, and molecules 
j2-4|. The study of confined atoms provides insight into various properties of quantum 
nanostructures like quantum dots or artificial atoms . For more detailed discussion on 
these applications we refer the reader to review articles |7H9|. 

-lecently numerous studies on helium atom confined in an impenetrable spherical cavity 



10N20| and in a penetrable spherical cavity 2l|, |22|] as well have been reported in the 



literature. Helium atom being the simplest many-electron system, the confined version of this 
atom provides a lucid way to study the effect of confinement on the electron correlation which 
arises due to coulomb interaction between the two electrons and pauli exclusion principle. 
Besides helium atom, few studies on some more confined many-electron atoms up to neon 



atom have also been reported in the literature [11|, |23|, |2J] . We note here that most of the 
studies on confined helium atom mainly considered the ground-state electronic properties. In 
contrast not many studies have been carried out on the electric response properties like dipole 
polarizability and hyperpolarizablity of this system. To the best of our knowledge results for 



the linear dipole polarizability of confined helium atom were reported only in Refs. |4|. 123 1 . 
However, no study devoted to the calculation of hyperpolarizability of confined helium atom 
exists in the literature. The main aim of this paper is to carry out calculation of not only 
linear dipole polarizability but also nonlinear hyperpolarizability of confined helium atom 
and study the evolution of these quantities with the size of the cavity (or strength of the 
confinement). We wish to point out here that the calculation of both linear polarizabilty 
(a) and third-order hyperpolarizability (7) of confined hydrogen atom were carried out in 
Ref. [25]. It was shown that both a and 7 strongly depend upon the radius of confinement 
and moreover, 7 changes sign and becomes negative under strong confinement. Thus it will 
be interesting to find out when 7 of confined helium atom undergo reversal of sign. 
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In this work we carry out calculation of electric response properties by employing density 



functional theory (DFT) based variation-perturbation approach |26|. To carry out these 
calculations we need to have ground-state densities of confined helium atom. This task 
has been accomplished by employing a variational approach involving minimization of the 
ground-state energy fuctional within the realm of DFT. A brief description of the theoretical 
methods employed for calculations of both ground-state densities and electrical response 
properties are presented in Section 2. The Section 3 is devoted to the discussion of results 
and paper is concluded in Section 4 

II. METHOD OF CALCULATIONS 



We begin this section with a brief description of the method for obtatining ground-state 
density of confined helium atom within the realm of DFT. The Kohn-Sham (KS) equation 
of DFT are obtained by minimizing the energy functional (in atomic units) j27| 



Eks[p\ = Ts[p] + J[p] + E,,[p] + j Ve,,{r)p{r)dv. (1) 

Here Ts[p] denote the kinetic energy functional of non- interacting particles and in terms of 
single-particle orbitals '?/'i(r) it is represented as 

N 

Ts[p] = Yl / ^;(r)(--V^)^.(r)dr. (2) 

i=i 

These orbitals yield density of interacting system via 

N 

i=l s 

In Eq. (II]) J[p] represent the classical part of the electron-electron repulsion, Ea;c[p] is the 
exchange-correlation functional and the last term corresponds to the contribution due to the 
external potential Vext{^)- For confined atom fext(r) has two parts namely: (1) the nuclear 
potential —Z/r (where Z is the nuclear charge) and (2) the confining potential Vconf due to 
an impentrable spherical box of radius Tc of the form 



Vconf {r) = < ' (4) 
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The minimization of Ei^s[p] with respect to single-particle orbital ipi{r) satisfying the con- 
dition given by Eq. ([3]) leads to so-called Kohn-Sham equation which is the workhorse of 
DPT. In this paper we carry out this minimization explicitly by using appropriate varia- 
tional forms for the single-particle orbitals. For this purpose we expand the single-particle 
orbital of a confined helium atom as 

= Yl ^iXi{r)fc{r) (5) 

i 

where {cj} are the variational parameters which are determined by minimization of the 
ground-state energy (Eq. ([1])), fc(r) is the cut-off function which takes care of the confine- 
ment boundary condition of density vanishing at r = rc, and the basis function is given 
by the product of a Slater-type orbital (STO) for the radial part and a spherical harmonic 
function Yim{0, 0) for the angular part as 

x{T) = RUr)Yim{9,(P). (6) 

The radial function Rri/(r) is given by 

Rni{r)=r^-'e-^'- (7) 



with n and ( representing orbital parameters which we choose from Ref. For cut-off 

function fJr) we choose both linear 



and the quadratic 



fcir) = ( 1 - - ) (8) 



fcir) =1-3 (9) 



forms and investigate their performance in yielding ground-state energies of confined helium 
atom. 

Having described the variational approach for obtaining ground-state density of confined 
helium atom we next briefly outline the method adopted in this paper to calculate static 
linear (a) and nonlinear (7) polarizabilities. The response properties mentioned above are 
calculated by employing variation-perturbation (VP) approach within DPT. In density based 
VP, energy to order (2n -|- 1) is determined by the perturbation expansion of the density 
correct to order n only. Purther, the even-order energy correction is minimum for 
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the exact induced density if expansion up to order n is known exactly. For details of 

the VP approach within DFT, we refer the reader to reference j26|. For our purpose here 
it is sufficient to note that a and 7 are calculated from the second-order /S.E^'^^ and the 
fourth-order AE^^'^ change in energies respectively, by employing relations 



7 = -24AE(^) (10) 



These energy changes are in turn obtained variationally by minimising 

E(^) = / v^'\r,)p('\r^)dr^ + \ f /^fw J '\n)p'^'\^2)dv,dv,^ (11) 
J 2 J 6p{ri)6p{r2) 

with respect to p^^\ and 

2 J 6p[vi)dp[Y2) 

+ \ i A i i / '\r,)p^'\r2)p('\r,)dr,dr2dr, 

2 J dp{ri}dp{r2)6p{r3} 



24 J dp{ri)dp{r2)dp{rs)dp{r4) 
X (iri(ir2(ir3(ir4. (12) 

with respect to p^'^\ Here v'^^)(r) is the applied (external) perturbation. F[p] is a universal 
functional of the density and it is given by the sum of the kinetic, Hartree and the exchange- 
correlation energies of the electrons. All the functional derivative in the equations above 
(Eqs. ffTT]) and (fT^ ) are evaluated at the ground-state density po- For aii atom placed in a 
static electric field £ along z-axis the variational ansatz for p^^^ and p^^^ are 

p(i)(r) = Ai{r)cos6po{r), 

p(2)(r) = [A2(r) + A3(r)cos2^]po(r) + Ap(°)(r) (13) 

where 

Ai{r) = air + hr^ + ar^ + ■■■ ,i = 1,2,3 ■■ ■ (14) 

with a.i, hi ■ ■ ■ being the variationals parameters. A is fixed for each set of parameters 
by the second-order normalization condition J p^'^\r)dr = 0. Notice that the first-order 
normalization condition J p^^\r)dv = is automatically satisfied by pW(r) in Eq.(fT3|). 
We have used five parameters for Ai and eight parameters each for A2 and A3. Adding 



more parameters does not affect the results significantly indicating their convergence. For 
evaluating the functional derivatives of the the exchange and correlation energies, we use 

for the exchange contribution and Gunnarsson-Lundquist (GL) 



the Dirac functional 



parametrization j29l] for the correlation energy within local-density approximation (LDA). 
In the next section we discuss results obtained by us using above-mentioned methods. 

III. RESULTS AND DISCUSSIONS 

We begin this section with the discussions of the results for ground-state energy of a 
confined helium atom obtained by us to assess their accuracies. In this connection we note 
that DFT based results for confined helium atom have already been reported in Ref. 
which were obtained by numerically integrating the KS equation with Dirichilet boundary 
condition 30| . In order to establish the accuracy of our variational results we compare them 
with those of Ref. jl^. First of all we note that we perform calcualtion with both linear 
and quadartic cut-off functions as given by Eqs. dH]) and respectively. We find that the 
ground-state energies for several values of confinement radius Vc obtained with quadratic 
cut-off function are close but slightly lower than the corresponding results obtained with 
linear cut-off functions. Therefore, in the following we report results only with quadratic 
cut-off functions. 

In table I we present the results for the energies for the ground-state ^S'(ls^) of a confined 
helium atom as a function of confinement radius Tc- In this table we present the results for 
the case of exact exchange (EXX) energy, which are obtained by substituting Ex = —Eh/ 2 
(exact for two-electron systems) and Ec = 0. This case corresponds to Hatree-Fock (HF) 



approximation and we compare our EXX results with those of Ref. We also present 

the results obtained with exchange-only (XO) with E^; = and exchange-correlation (XC) 
energies within LDA in second and third columns respectively. These results are compared 
with the corresponding numbers of Ref. which were obtained by numerically solving 
the Konn-Sham equation with Dirchilet boundary condition. In order to assess the accuracy 
of DFT based results we also display results obtained via correlated wavefunction based 



calculation with 7-parameter Hylleraas expansion 18[ in the last column of Table I. First 



we note that the results for the case of EXX energy obtained by us match very well up to 



4-th decimal place with the results of Ref. ll| for all values of r^. This establishes the 
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accuracy of the variational method employed by us. The XO-LDA results obtained by us 
are close but slightly higer than the corresponding EXX values as long as Vc < 1.0 a.u.. 
On the other hand, for > 1.5 a.u. we find that trend is just reverse. Moreover, XO- 
DA results are also quite close to the data available for the range Tc = 2 — 6 a.u. in Ref. 



17j . With the inclusion of correlation energy term within LDA the ground-state energies of 



confined helium atom reduce slightly as compared to the corresponding XO numbers. Our 



XC results match well with those of Ref. 



171] and the difference in the two results are mainly 



due to the use of two different XC functionals for the calculations. We note here that in 

n n 

Ref. [17| Pewrdew-Wang form of the correlation functional [311] along with Dirac form for 
the exchange energy functional has been employed whereas we employ GL parametrization 



for the correlation part 



29| . With the inclusion of correlation, the ground-state energy of 



a confined helium atom decreases relative to the corresponding XO-LDA valuses as long as 
energies remain positive. For confinement radii with negative ground-state energies inclusion 
of correlation leads to lowering of the ground-state energy. Similar trend is also observed 
with Perwdew-Wang XC functional. 

The comparison of DFT based results with the corresponding Hylleraas wavefunction 
based numbers clearly shows the EXX results are the closest to the latter. From these 
results we conclude that the contribution of correlation energy both in strong {vc < 1 a.u.) 
and weak (rc > 1 a.u.) confinemnets is not very significant to the total energy of a confined 
helium atom. Thus for confined helium atom it is possible to get sufficiently accurate results 
for the ground-state energy provided exchange part of the energy is accurately taken into 
account. 

Having established the accuracy of the DFT based results for the ground-state energy of a 
confined helium atom we now present the results for static polarizability a and second-order 
hyperpolarizabilty 7 of confined helium atom by employing the ground-state densities which 
are obtained via above-mentioned variational ground-state energy calculations. In Table 
11 we compile the results for both a and 7 corresponding to several values of Vc obtained 
with EXX, XO-LDA and XC-LDA ground-state densities. It can be seen from Table II 
that for all the three energy functionals a decreases monotonically with decrease in and 
appproaches zero for very small value of Tc. This trend is in agreement with the results 



of Ref. 



23|. We also note that with increase in Vc the values of a correctly tend to the 



respective free atom cases which are aExx = 1-33 a.u., axo = 1-77 a.u., and axe = 1-63 



7 



a.u. [32']. It is also interesting to note from Table II that for Tc < 1.5 a.u. the values of 
polarizability a obtained with three differnt densities are almost identical. Therefore, like 
ground-state energy of a confined helium atom its polarizability too does not have a strong 
dependence on the correlation energy specially in the strong confinment regime. Next we 
focus our attention on the results for the hyperpolarizability 7 of a confined helium atom, 
which are to our knowledge not reported earlier. The results for 7 also correctly converge 
to their respective free atom values with the increase in Vc- These values are 'Jexx = 36.2 
a.u., ■jxo = 114.2 a.u., and ■jxc = 88.15 a.u. [3^. Like polarizability, the values of 7 also 
show a monotonic decreasing trend, however, the decrease in the values of 7 are more rapid 
as comapred to the polarizablity. With the increase in the compression, that is decrease 
in Tc the values of hyperpolarizability 7 not only approach zero but also undergo a change 
in sign. The change in sign occurs at different value of Vc for three different ground-state 
densities employed in this paper. The change in the sign of hyperpolarizabilty of a strongly 
confined hydrogen atom has already been discussed in Ref. 251]. We note here that the 
chaiige in sign of 7 is akin to the term hypopolarizability which was coined by Coulson et 
al. 331]. The change in sign of 7 was later discussed by Langhoff et al. 3J] in connection 
with study of hyperpolarizability of high Z ions isoelectronic with Na and Mg series. They 
concluded that the compact electronic charge density resulting from the increasing value of 
Z within isoelectonic series is responsible for reversal of sign of 7. Similarly in a confined 
atom the charge density becomes highly compact with decreasing radius of confinement 
thereby yielding a negative hyperpolarizability. Finally, we note that unlike polarizabilty 
the results for 7 obtained with different densities vary for all values of Vc. This indicates 
that hyperpolarizability 7 depends crucially on the nature of ground-state density employed 
for the calculation. 



IV. CONCLUSION 

In this paper we have calculated the ground-state energies and the densities of a helium 
atom confined in an impentrable spherical box for various values of radius of confinement 
within DFT. We use three different exchange-correlation energies namely, exact exchange 
equivalent to HP level, exchange-only at LDA level, and correlation included with exchange 
at LDA level to obtain the ground-state properties. Using these ground-state densities we 
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perform calculations of static linear and nonlinear electric response propertise of confined 
helium atom and study their variation with the strength of the confinment. The ground-state 
energies and densities are obtained by variationally solving the KS equation of DFT with 
the variational form for the ground-state orbital expanded in terms of STOs. The results 
obtained by us are quite accurate and compare well with the already published data. We 
find that energies of a confined helium atom are not affected by the inclusion of correlation 
energy especially in the strong confiment regimes. With the increase in the value of radius 
of confinement the differences between the results for energy obtained by employing three 
different densities increase indicating the importance of correlation energy for determining 
the energy of a free helium atom. The linear polarizability a of a confined helium atom 
decreases with the decrease in the value of confinement radius approacing the value of zero 
in the strong confinement regime. Results of our calculations clearly demonstrate that for 
Tc < 1.5 a.u. the values of a do not depend much on the correlation energy. On the other 
hand, value of hyperpolarizability 7 shows a decreasing trend with decreasing Tc which 
approaches the value of zero rapidly and then changes sign on further decreasing Tc. This 
change in sign of 7 is attributed to charge density becoming highly compact in the case of 
strong confinement. 
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TABLE I: Ground-state energies of confined helium atom as a function of confinment radius 
(in atomic units) for the case of exact exchange (EXX), exchange-only within LDA (XO-LDA) 
and exchange-correlation within LDA (XC-LDA). The numbers in the parenthesis are taken from 
already published data. The results in the last column are taken from 18 1 



rc 


EXX 


XO-LDA 


XC-LDA 


Correlated- Hylleraas 




(Ref.fll]) 


(Ref. [17]) 


(Ref. m 


(Ref. [18]) 


0.5 


22.79096 


23.32202 


23.099 


22.7419 




(22.79095) 


- 


- 




0.6 


13.36683 


13.81792 


13.605 


13.3187 




(13.36682) 


- 


- 




0.7 


7.97302 


8.36716 


8.164 


7.9258 




(7.97302) 


- 


- 




0.8 


4.65736 


5.00895 


4.8133 


4.6110 




(4.65737) 


- 


- 




0.9 


2.50942 


2.82805 


2.369 


2.4638 




(2.50944) 


- 


- 




1.0 


1.06120 


1.35362 


1.17040 


1.063 




(1.01624) 


- 


- 




1.5 


-1.86422 


-1.64897 


-1.81185 


-1.9066 




(-1.86422) 


- 


- 




2.0 


-2.56257 


-2.39363 


-2.53480 


-2.6035 




(-2.56253) 


(-2.38363) 


(-2.50589) 




3.0 


-2.83078 


-2.68201 


-2.82256 


-2.8715 




(-2.83083) 


(-2.68210) 


(-2.79608) 




4.0 


-2.85854 


-2.71807 


-2.85552 


-2.8994 




(-2.85852) 


(-2.71813) 


(-2.82970) 




5.0 


-2.86138 


-2.72288 


-2.85856 


-2.9026 




(-2.86134) 


(-2.72290) 


(-2.83387) 




6.0 


-2.86162 


-2.72346 


-2.86000 


-2.9032 




(-2.86151) 


(-2.72354) 


(-2.83439) 
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TABLE II: Static polarizability a and hypcrpolarizability 7 of a confined helium atom for various 



values of confinemnt radius Vc- AH numbers are in atomic unit 



1'c 


EXX 


XO-LDA 


XC-LDA 


(\ 




n 




n 




1.0 




-± .UO A lU 




1 lRvin~4 




S 78 V 1 n"~5 
-0. /OA iU 


1.1 


u.uoy 


1 70N^in^4 
-1. / Z X iU 


U.UDU 


-Z.U4 X lU 


U.UDU 


-i.oDX lU 


1.2 


u.u ( y 


9 1 c 1 n— 4 
-Z.lOX lU 


n nsn 

U.UoU 


-z.y i X iu 


U.UoU 


-l.^DX lU 


1.3 


n 1 no 


1 0/1 s/ 1 n~4 


U.iUo 


-O.IZX iU 


U.iUo 


ns V 1 n~6 

-Z.Uo X iU 


1.4 


u.izy 


O.OftX iU 


U.iOi 


7 0/1 V 1 n~5 


U.iOi 


c; oos/ 1 

O.ZZ X lU 


1.5 


U. iOU 


1.00 A lU 


n 1 RO 


7 7/L V 1 n~4 


n 1 RO 

U. iOZ 


l.c5 / X lU 


2.0 


0.360 


0.0603 


0.374 


0.051 


0.372 


0.0634 


2.5 


0.606 


0.487 


0.647 


0.469 


0.641 


0.530 


3.0 


0.839 


1.900 


0.934 


2.158 


0.916 


2.324 


4.0 


1.161 


10.358 


1.392 


14.862 


1.336 


14.674 


5.0 


1.286 


23.381 


1.641 


45.029 


1.544 


40.552 


6.0 


1.322 


35.766 


1.763 


93.056 


1.628 


75.265 
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